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Abstract 

The proprotein convertase site-1 protease (SIP) converts latent ER-membrane bound transcription factors SREBPs and ATF6 
to tlieir active forms. SREBPs are involved in cholesterol and fatty acid homeostasis whereas ATF6 is involved in unfolded 
protein response pathways (UPR). Cartilage-specific ablation of SIP in mice (S1P^''°) results in abnormal cartilage devoid of 
type II collagen protein (Col II). S1P'^''° mice also lack endochondral bone development. To analyze S1P'^''° cartilage we 
performed double-labeled immunofluorescence studies for matrix proteins that demonstrated that type IIB procollagen is 
trapped inside the ER in S1P'^''° chondrocytes. This retention is specific to type IIB procollagen; other cartilage proteins such 
as type IIA procollagen, cartilage oligomeric matrix protein (COMP) and aggrecan are not affected. The SI P'^''° cartilage thus 
exhibits COMP-, aggrecan-, and type IIA procollagen-derived matrices but is characterized by the absence of a type IIB 
procollagen-derived matrix. To understand the molecular reason behind S1P'^''° phenotypes we performed genome-wide 
transcriptional profiling of cartilage isolated from SIP'^''" and wild type littermates. While the UPR pathways are unaffected, 
the SREBPs-directed cholesterol and fatty acid pathways are significantly down-regulated in S1P'^''° chondrocytes, with 
maximal down-regulation of the stearoyl-CoA desaturase-1 (Scdl) gene. However, mouse models that lack Scdl or exhibit 
reduction in lipid homeostasis do not suffer from the ER retention of Col II or lack endochondral bone. These studies 
indicate an indispensable role for SIP in type IIB procollagen trafficking from the ER. This role appears not to be related to 
lipid pathways or other current known functions of SIP and is likely dependent on additional, yet unknown, SIP substrates 
in chondrocytes. 
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introduction 

Site-1 protease (SIP; also called the membrane-bound tran- 
scription factor protease, site-1) is a proprotein convertase that 
converts latent, endoplasmic reticulum (ER) membrane-bound 
transcription factors into their free and active form. Two 
developmental pathways regulated by SIP that have been studied 
extensively include cholesterol and fatty acid homeostasis and the 
unfolded protein response [1]. During cholesterol and fatty acid 
homeostasis, S 1 P plays a fundamental role in the processing of the 
transcription factors sterol regulatory element binding proteins 
(SREBP-la, -Ic, and -2) [2]. During unfolded protein response 
(UPR), SIP plays a critical role in the processing of activating 
transcription factor 6 (ATF6) [3], old astrocyte specifically induced 
substance (OASIS) [4], and the cAMP-responsive element binding 
protein H (CREBH) [5]. All these pathways are fundamental in 



maintaining cellular homeostasis and therefore SIP plays major 
and critical roles in fundamental developmental pathways. 
Mutational inactivation of SIP in zebrafish (gonzo) results in 
abnormal lipid metabolism and skeletal formation [6]. 

Mammalian skeletal development begins early during embryo- 
genesis. It involves endochondral bone formation, a process by 
which most long and axial bones of the body are formed [7]. 
During this process, mesenchymal cells aggregate to form 
condensations and cells in the center of these condensations 
differentiate to form chondrocytes [8]. Chondrocytes secrete a 
specialized extracellular matrix (ECM) called cartilage that forms 
the anlagen for future skeleton, whereas the surrounding 
undifferentiated cells form the perichondrium [9]. The innermost 
cells within the chondrocytes differentiate to form hypertrophic 
chondrocytes, with cells in the inner layer of the perichondrium 
differentiating into bone-forming osteoblasts resulting in a bone 
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collar around the cartilaginous core [10-12]. Fully difTerentiated 
hypertrophic cells secrete a distinct ECM which gradually 
becomes calcified before they die through apoptosis. The 
mineralization of the ECM allows for vascular invasion from the 
bone collar and the entry of osteoclasts and osteoblasts that 
degrade the mineralized cartilage to deposit endochondral bone 
[13]. Deposition of a type I collagen-rich bone matrix results in 
two opposing growth plates allowing for bone growth in both 
directions [14]. In contrast to endochondral bone formation, much 
of the craniofacial skeletal bones are formed by intramembranous 
ossification where the mesenchymal cells diflFerentiate directly into 
osteoblasts without requiring a cartilage intermediate. 

It follows from this developmental paradigm that factors that 
affect cartilage development also impact endochondral bone 
development. The link of SIP to cartilage development was 
shown by the study of zebrafish gonzo mutant [6] . In a previous 
study we showed that SIP is required for proper cartilage matrix 
development in mice [15]. By creating cartilage-specific SIP 
knockout mice (SlF*°) we demonstrated that these mice do not 
form any endochondral bone. S1P*° mice also exhibited poor 
cartilage development with most of the type II coUagen protein 
(Col II) trapped inside the cell, resulting in a drastic reduction of 
Col II in the cartilage. Ultrastructural analysis of the cartilage 
showed engorged and fragmented ER. 

In the current study we investigated the nature of Col II 
entrapment and the mechanistic reasons behind SIF*" pheno- 
types. Lack of SIP would result in lack of activation of SREBPs, 
ATF6, OASIS, and CREBH. Therefore, lack of SIP activity in 
chondrocytes would be expected to affect both the SREBPs- 
directed cholesterol and fatty acid homoeostasis and the UPR 
pathways. In order to understand how lack of SIP affects the 
downstream pathways, transcriptional profiling in chondrocytes 
was performed by genome-wide expression analyses with RNA 
extracted from the cartilage of S1P^*° and wild type (WT) 
Httermates. Our studies show that the SREBPs-dependent 
cholesterol and fatty acid biosynthetic pathways are down- 
regulated in SIP^*" chondrocytes. In contrast, UPR pathways 
remain unaffected. Furthermore, lack of SIP in cartilage results 
specifically in the ER retention of type IIB procollagen (pro-Col 
IIB). These data suggest that SIP has an indispensable function in 
pro-Col IIB traflicking from the ER to the cartilage matrix. 
However, our in depth mechanistic analyses indicate that this 
activity is not related to current known functions of SIP. 
Additional, yet unidentified, SIP substrates presumably modulate 
Col II trafficking from the ER to the cartilage ECM. 

Materials and Methods 

Ethics Statement 

All mouse procedures were performed in accordance with 
National Institutes of Health's Guide for the Care and Use of 
Laboratory Animals using vertebrate animals /ethics protocols 
reviewed and approved by the Animal Studies Committee at 
Washington University School of Medicine. 

Double-labeled immunofluorescence studies 

Double-labeled immunofluorescence to analyze retention of 
matrix proteins in the ER were performed on 5- |xm formalin 
fixed, paraffin embedded sections from El 6.5 embryonic tissues 
processed as described [15]. Following rehydration through 
graded ethanol solutions, the sections were washed in phosphate 
buffered saline, pH 7.4 (PBS) for 10 min. For antigen retrieval the 
slides were treated with Proteinase K (10 |ig/ mL in 1 0 mM Tris- 
HCl pH 7.4) for 30 min at 37°C followed by blocking for 1 h in 



either 10% donkey serum or 10% goat serum. The two primary 
antibodies (Table 1) were added together in 1.5% donkey serum 
(or goat serum) and incubated in a humidified chamber for 1 h. 
The slides were washed three times in PBS for 15 min and rinsed 
once in distilled water. The two secondary antibodies were then 
added together and incubated for 1 h, washed thoroughly and 
mounted with Vectashield DAPI (Vector Laboratories, Inc., CA). 
Antibody concentrations were as follows: antibodies to Col II triple 
helical domain (THD; IIF antibody [16]), Col ILA. (designed 
against exon 2; IIA antibody) [15,16], aggrecan (rabbit antisera 
against rat aggrecan, a gift of Dr. Kurt Doege), and cartilage 
oligomeric matrix protein (COMP) [17] were used at 1:100. 
Affinity-purified IIBN antibody (designed against type II coUagen 
exon 1-3 peptide junction and specific to type IIB procollagen) 
[18] was used at 5 ng/mL, while goat anti-human calnexin 
(Santacruz Biotechnology) and anti-XBP-1 (Biolegend, CA, USA) 
antibodies were used at 1:50. The secondary donkey anti-goat 
Alexa 488 (Invitrogen) and donkey anti-rat Alexa 594 (Invitrogen) 
antibodies were all used at 1:250. Images were captured using a 
60X, 1.4 NA oil immersion objective mounted on an Eclipse E800 
microscope (Nikon) and Qlmaging Retiga 2000R Fast 1394 
camera and deconvolved. For deconvolution imaging, Meta- 
Morph software (Molecular Devices) was used to control the Z- 
motor device (Prior Scientific), to capture and deconvolve images, 
and to compile them to give a final image. All images shown are 
for mature columnar proliferative chondrocytes in the WT and for 
a corresponding region in SlI*^*" (Fig. SIA). 

Genome-wide expression profiling 

To profile gene expression in SlP*" cartilage and to contrast it 
with WT, we harvested the chondroepiphyseal cartilage from the 
tips of long bones in foreUmbs/hindlimbs (Fig. SIB) of freshly 
harvested embryonic (E) 16.5 SIP"*" (SlP'^^;Col2-Cre) embryos 
and its WT (SIP^, SIF''"', or SlP^'''";Col2-Cre) Httermates, as the 
mutant mice die during or very soon after birth. At El 6.5 the 
mutant is easily distinguished without genotyping, plus identifying 
and separating the bones from the musculature is easier than at 
later time points in the mutant. After harvesting the bones, excess 
muscle tissue and blood vessels were removed by rolling on 
Whatman filter paper, and the chondroepiphyseal cartilage from 
the ends of the bones were broken off (Fig. SIB) and stored in 
Trizol (Invitrogen) at — 80°C for RNA extractions later. In some 
cases the chondroepiphyseal cartilage was not separated and the 
entire long bones kept intact and stored in Trizol for RNA 
extraction. The calvarium was also surgically removed in some 
cases and stored in Trizol at — 80°C and RNA extracted as for the 
chondroepiphysis. Therefore three kinds of tissue were harvested: 



Table 1. Summary of antibodies used for 
immunohlstochemical analyses. 





Antibody 


Protein recognized 


IIF 


Type II collagen triple helical domain (THD) 


IIA 


Type IIA procollagen (Col IIA} 


IIBN 


Type IIB procollagen (pro-Col IIB) 


Anti-calnexin 


Calnexin 


Anti-XBP-1 


XBP-1 


Anti-aggrecan 


Aggrecan 


Anti-COMP 


COMP 
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chondroepiphyseal cartilage from the tips of long bones, entire 
long bones (includes both the bone and epiphyseal cartilage), and 
calvaria. 

To extract RNA, the harvested tissues were first homogenized 
in Trizol using a polytron and RNA extracted according to the 
manufacturer's recommended protocol and column purified by 
RNeasy mini kit (Qiagen). Genome-wide expression profiling was 
performed using lUumina MouseWG-6 vl.l Expression BeadChip 
which profiles 45,200 mouse transcripts. Statistically significant 
difierentially expressed genes (false discovery rate (FDR) of 5%) in 
SIP'*" chondrocytes as compared to WT were identified by 2-way 
ANOVA (analysis of variation) analysis using Partek Genomics 
Suite (St. Louis, MO). Only genes that were two-fold or more 
differentially regulated with significant j!>-values {p<0.05) were 
carried forward for characterization. DifiTerential gene expression 
data derived from microarray analysis was used to identify 
functional and molecular networks through the use of MetaCore 
GeneGo network building tools (GeneGo Inc., Carlsbad, CA). 
Again, only genes that were two-fold or more differentially 
regulated were used to generate a network score of negative log of 
the p value to determine biological pathways disrupted in SIP'*" 
chondrocytes. The raw microarray data have been deposited in 
the Gene Expression Omnibus (GEO). The accession number is 
GSE55577 (online at http://www.ncbi.nlm.nih.gov/projects/ 
geo). 

Quantitative real-time PCR (qPCR) analysis 

Quantitative real time PCR (qPCR) was performed for 
confirmation of expression profiling for genes involved in 
cholesterol and fatty acid homeostasis. For qPCR analysis, the 
RNA was harvested and pooled from the chondroepiphyseal 
cartilage of forelimbs and hindlimbs of five E16.5 S1P^*° and five 
E16.5 WT (three SIP'"'^, one SIP'^^ and one Sl/'''';Col2-Cre 
mice) littermates. The RNA was extracted with Trizol and column 
purified as described above and used for qPCR analysis using 
SYBR Green primer sets, 2 x SYBR Green mix (Life Technol- 
ogies/Applied Biosystems) and standard protocols. SYBR Green 
primer sets for the murine cyclophilin, insulin-induced gene-1 
(Insig-1), 3-hydroxy-3-methylglutaryl (HMG)-CoA reductase 
(Hmgcr/Hmdh), HMG CoA synthase (Hmgcs). farncsyl diphos- 
phate synthase (Fdps), low density lipoprotein receptor (Ldlr), 
stearoyl-CoA desaturase-1 (Scdl), and Scd2 genes have been 
reported previously [19,20]. Additional primer sets for the murine 
glyceraldehyde 3-phosphate dehydrogenase (Gapdh), fatty acid 
desaturase-2 (Fads2), sterol-C4-methyl oxidase-like (Sc4mol), star- 
related lipid transfer (START) domain containing 4 (Stard4), 
proteoglycan 4 (Prg4), loricrin (Lor), and fibroblast growth factor 
receptor 1 oncogene partner 2 (Fgfrlop2) genes were taken from 
the MGH/Harvard Medical School primer bank (http://pga. 
mgh.harvard.edu/primerbank/) and are given in Table SI. The 
relative amount of all mRNAs was calculated using the 
comparative C, method with cyclophilin or Gapdh as the invariant 
control. 

Analysis of ER stress response in Sip'^'^° chondrocytes 

RNA extracted from the chondroepiphysis of El 6.5 WT and 
SIP'*" embryos was used to further validate the unaffected UPR 
discerned from the microarray analysis. For X-box binding 
protein- 1 (XBP-1) mRNA splicing assay, RNA was converted to 
cDNA using routine protocols and splicing analyzed by PCR or by 
restriction digestion. For PCR analysis, the cDNA was amplified 
using XBP-1 primers as described [21] and the PCR products 
directly analyzed in a 4.8% polyacrylamide gel to separate the 
products derived from the spliced and un-spliced XBP-1 mRNA 



based on size. Analysis by restriction digestion was performed as 
described [22]. Briefly, the cDNA was first amplified using the 
XBP-1 primers 3S and 12AS followed by digestion of the 
amplified PCR product with Psll to selectively cleave the cDNA 
derived from un-spliced XBP-1 mRNA, and analyzed in a 2% 
agarose gel. 

For a global analysis of ER stress response in S1P"*° 
chondrocytes the RT^ Profiler PCR Array System for the murine 
UPR (Cat. No. PAMM-089C-2, SABiosciences, Qiagen), which 
profiles the expression of 84 key genes involved in recognizing and 
responding to unfolded protein accumulation in the ER, was used 
(a complete list of proteins profiled can be found at http:// 
sabiosciences.com/rt_pcr_product/HTML/PAMM-089A.html) 
(also listed in Table S2). The RNA was converted to cDNA using 
the RT^ First Strand kit (Cat. No. C-03) and qPCR reactions set 
up using the RT^ SYBR Green/ROX qPCR master mix (Cat. 
No. PA-012) in a 96-well array plate following instructions 
recommended by the manufacturer. A total of four qPCR arrays 
were used, two for WT and two for SIP'*" RNA for biological 
duplicates. Each WT or SlF*" RNA was pooled from the 
chondroepiphyseal cartilage of two E16.5 WT or two E16.5 SIP'*" 
embryos, respectively. Thus a total of four WT and four SIP'*" 
El 6.5 embryos were profiled. Analysis of the data was done using 
the Qiagen-SABiosciences' web-based PCR data analysis software. 

Miscellaneous 

mice were generated by mating S 1 P'^ mice witii SIP^''"'; 
Col2-cre mice as described [15]. Scdl — /— mice were generated 
by mating Scdl+f — heterozygotes and identified by genotyping 
[23]. Scdl+/— heterozygote parental strains and the SIP''^ and 
SlP'''*';Col2-Cre strains were all fed ad libitum with standard 
laboratory chow. Analysis of endochondral bone formation in 
Scdl—/— mice by Safranin O, Fast green, and hematoxylin 
staining and Col II deposition by double-labeled immunofluores- 
cence were done as for SIP'*" mice at El 5.5. In situ hybridization 
analyses w("re done on 5- |j.m paraffin-embedded sections as 
described previously using "^''S-labeled riboprobes [15,24]. Ribop- 
robes for the murine Scdl, Scd2, Ldlr, Fads2, and Fdps genes 
were derived from their full-length cDNA clones purchased from 
commercial sources. Murine Scdl immunohistochemistry (IHC) 
was performed on 5- |J,m paraffm-embedded sections using anti- 
Scdl antibody (Santacruz Biotechnology) and HRP-conjugated 
anti-goat IgG and DAB. Hematoxylin was used as the countc-r 
stain. In situ hybridization images were viewed with BX51 
(Olympus) microscope and images captured with a digital camera 
(DP70; Olympus) using DP controller software (Olympus). Images 
of hybridization signals were artificially colored red and superim- 
posed on toluidine blue-counterstained images using Photoshop 
(Adobe). 

Results 

Retention of cartilage matrix proteins in the ER 

Chondrocytes by nature are secretory cells and secrete large 
amounts of ECM proteins needed to create the cartilage matrix. 
Ultrastructural analysis of cartilage from SIP'*" mice showed an 
enlarged ER with the retention of a crystalline material [15]. 
Therefore we first examined the composition of proteins retained 
within the ER. Using double-labeled immunofluorescence, an 
antibody to calnexin (an ER membrane protein) and an antibody 
(IIF; see Table 1) against the Col II triple helical domain (THD), 
we analyzed whether Col II in the mutant cartilage was trapped in 
tiie ER (Fig. 1). 
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Calnexin Col II Calnexin + Col II 




XBP-1 Col II XBP-1 + Col II 




Figure 1. The ER retention of Col II in SIP chondrocytes. (A-F) Double-labeled immunofluorescence analyses for Col II THD (IIP antibody) 
and calnexin in El 6.5 femurs in WT (A-C) and SI P"^*" (D-F). Signals from calnexin (A, D) are shown in green (white arrows) and those from Col II THD 
(B, E) are shown in red (white arrowheads) with composite signals shown in C and F. (G-L) Double-labeled immunofluorescence analyses for Col II 
THD and XBP-1 in El 6.5 femurs in WT (G-l) and SIP"^*" (J-L). Individual signals from XBP-1 (G, J) are shown in green (white arrows) and those from Col 
II THD (H, K) are shown in red (white arrowheads), with composite signals shown in I and L. All images shown are for mature columnar proliferative 
chondrocytes in the WT and a corresponding region in S1P'^''° (see Fig. SI A). Bar (all panels): 10 |im. 
doi:1 0.1 371 /journal.pone.01 05674.g001 



In WT, the Col II THD in the cartilage (shown as red 
immunofluorescence, white arrowheads. Fig. IB, C) is distinct 
from calnexin (shown as green immunofluorescence, Fig. 1 A, C) of 
the ER membrane surrounding the blue DAPI-stained nuclei of 
the chondrocytes (Fig. lA-C). Notice the strong presence of a Col 



Il-derived matrix in the cartilage of the WT (Fig IB, C). In S1P'*° 
however, signaling from Col II THD (red) is visible primarily from 
inside the cells (white arrowheads. Fig. IE) concomitant with die 
drastic drop of Col ITderived signals from the SIP'*" cartilage. 
The chondrocytes appear distended from the entrapped Col II 
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protein. In the composite image (Fig. IF), Col II is visible as 
yellowish-orange immunofluorescence signal surrounding the 
blue-DAPI-stained nuclei, due to its colocalization with calnexin 
(white arrows, Fig. ID, F) in the ER. To confirm the ER- 
entrapment of Col II in the mutant, we analyzed whether Col II 
colocalized with XBP-1 [25], a known ER luminal protein 
(Fig. IG-L) in Sir*" cartUage. As shown in Fig. IJ-L, Col II 
colocalized with XBP- 1 in the mutant chondrocytes confirming its 
entrapment in the ER. Even though results are shown for the 
proliferative zone. Col II retention is seen throughout all zones of 
the growth plate and is not limited by cell size, morpliolog}" or 
location in the growth plate. These data indicate that the matrix in 
gjpc*o (.^j.jji^gg Js defective and devoid of Col II due to disruption 
of Col II trafficking from ER to the cartilage. 

During synthesis. Col II is made as a procollagen of three 
identical pro-ot-chains consisting of the amino (N) terminal 
propeptide, the THD domain, and the carboxy (C) terminal 
propeptide. The Col II protein incorporated into the mature 
cartilage is the THD domain with the N- and C-terminal 
propeptides removed. To understand the nature of the trapped 
Col II, we performed double-labeled immunofluorescence with the 
IIBN antibody [18], which detects the unique exon 1-3 peptide 
junction sequence specific to type IIB procollagen (pro-Col IIB) 
produced by chondrocytes (in contrast to the type IIA procollagen 
made by chondroprogenitor cells), in combination with anti- 
calnexin (Fig. 2A, B). In WT, signals from pro-Col IIB are 
undetectable presumably due to a rapid processing of the 
procollagen (Fig. 2 A) [18]. In the composite picture for the WT 
(pro-Col IIB + calnexin; Fig. 2A) only signals from calnexin are 
seen showing the presence of the ER compartments. In contrast in 
S1P^*°, significant signals from pro-Col IIB are detected that 
overlaps with signals from calnexin inside the cell (Fig. 2B). 

When we repeated the analysis for pro-Col IIB in combination 
with IIF antibody (Fig. 2C-F), only the Col II THD is seen in the 
WT cartilage matrix indicating the presence of the mature 
processed Col II protein (Fig. 2C, E). However, hardly any Col II 
THD is detected in the SlF*" cartilage (Fig. 2D, F; also see 
Fig. IE). Signals from Col II THD is seen primarily trapped inside 
the cells (Fig. 2D) that overlaps with intracellular signals from pro- 
Col IIB (Fig. 2F). These observations indicate that in the mutant, 
the pro-Col IIB isoform of Col II is trapped in the ER. In order to 
analyze if this retention is specific to pro-Col IIB, we repeated the 
analysis using the IIA antibody specific to type IIA procollagen 
(Col'lIA). In the WT (Fig. 2G~T), Col IIA is detected only outside 
the cell in the cartilage ECM where it forms an organized lattice 
network (in green, Fig. 2H). Importandy, an organized lattice 
network specific to Col IIA is also seen in SIP'*" cartUage. The IIF 
antibody as seen before did not identify any Col II network in the 
mutant (compare Fig. 21 to 2L). These observations indicate that 
the absence of SIP activit)^ affects only pro-Col IIB trafficking 
from the ER, and not Col IIA. 

Next, we analyzed if trafficking of other major matrix proteins 
such as COMP and aggrecan is affected in SIP"*" (Fig. 3). In WT 
chondrocytes, COMP (in red) is visible in the cartilage (white 
arrowheads, Fig. 3B, G) and its localization is distinct from 
calnexin (in green, white arrows. Fig. 3A,C) which surrounds the 
DAPI-stained blue nuclei. The distribution of COMP in the WT 
cartilage matrix is similar to that of Col II (compare Fig. 3B to Fig. 
IB or 2C, E). Therefore when we labeled both Col II THD (red) 
and COMP (green) in the WT (Fig 3G), the staining patterns for 
Col II and COMP overlapped in the matrix; the cartilage is seen 
as yellowish-green in most areas with very little Col II or COMP 
detected intracellularly. Interestingly, the SIP^*" cartilage also has 
a COMP-derived matrix (Fig. 3E, F). But COMP distribution 



does not resemble that seen in the WT cartilage and appears 
disordered and filamentous. In fact, COMP distribution in SIP'*" 
cartilage approximates the distribution of Col IIA (compare 
Fig. 3E to 2K). This is probably due to the absence of the Col II 
THD domain in the SIR*" cartilage. Labeling SIR*" cartilage for 
both Col II (red) and COMP (green) (Fig. 3H) showed that COMP 
and Col II occupy distinct domains. While most of Col II is 
detected intracellularly, the majority of COMP is seen as a 
network in the cartilage. While some COMP is also detected 
trapped in the ER in SlR*° chondrocytes (white arrows with 
asterisk. Fig. 3E, H), these are usually very small aggregates and 
are usually surrounded by the larger Col II aggregates (Fig. 3H) 
suggesting retention due to a physical interaction with Col II (see 
DISCUSSION). Moreover, this retention of COMP is insignifi- 
cant when compared to pro-Col IIB retention and has no 
consequence in the formation of a COMP matrix in the mutant 
cartilage. Co-localization studies with aggrecan (Age) show that 
like COMP, SIP ablation does not affect its secretion (Fig. 3J). 
Thus our studies demonstrate that pro-Col IIB is the primary 
cartilage matrix protein trapped in the ER of SIP^*" chondrocytes. 

Genome-wide expression profiling 

The primary, well studied role of SIP is its transcription factor 
processing activity. Lack of transcription factor processing by SIP 
would affect both UPR and cholesterol/fatty acid homeostasis, 
which could be causal to SIR*" phenotypes. To understand the 
mechanisms behind this phenotype we profiled gene expression in 
mutant cartilage by microarray and contrasted it to WT cartilage 
to identify the biological processes disrupted in SIR*" chondro- 
cytes. Statistically significant genes were identified by 2-way 
ANOVA analysis coupled with FDR of 5% (only 5% genes 
identified in the analysis could be false positive). Only genes that 
were 2-fold or more differentially regulated in the mutant were 
considered. Our primary aim was to understand the difference in 
gene expression in the chondroepiphyseal cartilage which at 
embryonic stages is the bulk of the cartilage in the skeletal 
elements. While SlR*° does not have any endochondral bone, it 
does have an enhanced cortical bone. Therefore entire long bones 
were also harvested and compared in this study. The lUumina 
Mouse WG-6 vl.l Expression BeadChip used in this study was 
chosen as it profiles more than 45,200 mouse transcripts including 
probes derived from the Mouse Exonic Evidence Based Oligonu- 
cleotide set and some less-characterized genes/probe sets derived 
from exemplar protein coding sequences from RIKEN FAN- 
TOM2. 

A principal components analysis (PCA) scatter plot representing 
only 40.8% of the whole genome information present in the 
microarray hybridization signals demonstrated that very httie 
difference in gene expression is seen between the WT and SIP^*" 
calvaria (Figure SIC). This suggested that skeletal structures 
arising from intramembranous ossification were likely not as 
affected by SIP ablation presumably because this developmented 
pathway does not need a cartilage intermediate. The PCA analysis 
also revealed that the largest difference in gene expression profile 
between WT and S1P^*° lay in the cartilage. A 2-way ANOVA 
analysis coupled with a minimum of +2 to -2 differential fold 
regulation generated a stringent list of genes differentially 
regulated in SIP'*" chondrocytes. Com|3arison of the gene 
expression profile between WT and SIP" " long bones showed 
similar data (not shown) presumably due to the presence of RNA 
from the cartilage, indicating that the difference lay primarily in 
the cartilage. Incorporation of calvaria data set in these analyses 
did not make any difference to the differentially expressed genes 
identified. Tables S3 and S4 list genes that are significandy 
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Figure 2. Specific retention of pro-Col IIB in SIP cKiondrocytes. Double-labeled immunofluorescence analyses were done for pro-Col IIB 
(red, IIBN antibody) and calnexin (green) in El 6.5 femurs in WT (A) and SIP'^'"' (B) and only the composite for this analysis is shown. In a separate 
analysis (symbolized by the dashed line between panels A/B and C/D), double-labeled immunofluorescence was done for Col II THD (red) (C-F) and 
pro-Col IIB (green) (E, F) in El 6.5 femurs in WT (C, E) and SIP'^'"' (D, F) with composite signals shown in E and F. Panels G-L show double-labeled 
immunofluorescence analyses for Col 11 THD (red) (G, J) and Col IIA (green; IIA antibody) (H, K) in El 6.5 femurs in WT (G-l) and SIP"*" (J-L). 
Composite signals are shown in panels I and L. All images shown are for mature columnar proliferative chondrocytes in the WT and a corresponding 
region in SIP^'"'. Bar (all panels): 10 [im. 
doi:1 0.1 371 /journal.pone.01 05674.g002 



differentially regulated (up- or down-regulated, respectively) in 
SIP^*" chondrocytes. Figure 4A shows a volcano plot of the 
ANOVA analysis representing genes that are significantly differen- 



tially regulated (as per their ^ values) and exhibit fold changes of 2 or 
more. Figure 4B shows the range of differential regulation in SIP^*" 
chondrocytes and the identity of some critical genes. 
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Figure 3. COMP and aggrecan are not intracellularly retained in SIP*^*" chondrocytes. Double-labeled immunofluorescence analyses for 
COMP and calnexin in El 6.5 femurs in WT (A-C) and SI P"^*" (D-F). Individual signals from calnexin (A, D) are shown in green (white arrows) and those 
from COMP (B, E) are shown in red (white arrowheads), with composite signals in C, F. The relative distribution of Col II THD (red) and COMP (green) 
proteins for WT (G) and SI 9^''° (H) are shown as composites from double-labeled immunolabeling analyses in El 6.5 femurs. Occasional entrapment of 
COMP is indicated by arrows marked with an asterisk. Composite signals showing the localization of aggrecan (Age) (red) relative to calnexin (green) 
in El 6.5 WT (I) or SIP*^*" (J) femurs. All images shown are for mature columnar proliferative chondrocytes in the WT and a corresponding region in 
SIP^*". Bar (all panels): 10 fim. 
doi:1 0.1 371/journal.pone.01 05674.g003 

Analyses of difTerentially expressed genes to identify functional 
biological pathways disrupted in the cartilage due to SIP ablation 
were done by GeneGo (MetaCore). Surprisingly, this analysis 
showed that the top 10 processes upregulated in SIP'*" belonged 



primarily to UPR pathways (Fig.5A). This suggested that ER stress 
response pathways were unaffected in SIP^*" cartilage. Further- 
more, when a functional molecular network was created with these 
up-regulated genes with the addition of ATF6a to the network 
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Figure 4. (A) A volcanic plot of the 145 differentially regulated genes in the S1P'^''° cartilage identified by 2-way ANOVA analysis. Only genes that 
show at least a 2-fold differential regulation with a p<0.05 (FDR <0.05) are shown, with the color red and its gradations identifying up-regulated 
genes and the color blue and its gradations identifying down-regulated genes. (B) Schematic representation of critical genes up- or down-regulated 
in SIP"^*" chondrocytes as compared to WT showing the range of differential expression in the mutant (un-abbreviated gene names are listed in 
Tables S3 and S4). 

doi:10.1371/journal.pone.0105674.g004 
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Figure 5. (A) A schematic representation of the top ten biological pathways up-regulated in SIP chondrocytes identified by GeneGo (MetaCore) 
based on their p value. The majority of up-regulated pathways belong to UPR in response to ER stress. (B) A molecular network generated from 
upregulated genes with the addition of ATF6o( to the network. Green lines represent activation and red lines represent inhibition. 
doi:1 0.1 371 /journal.pone.01 05674.g005 

(Fig. 5B), a number of genes were found to be connected to SIP^*" chondrocytes. This suggested tliat ATF6 activity was 
ATF6a, directly or indirectly, and also to each other, acting in unlikely to be disrupted in S1P^*°. When the significandy 
concert to respond to ER stress due to pro-Col IIB entrapment in differentially down-regulated genes were analyzed similarly by 
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Figure 6. (A) A schematic representation of the top ten biological pathways down-regulated in SIP'^'"' chondrocytes as identified by GeneGo based 
on their p value. The majority of down-regulated pathways identified belong to cholesterol and lipid biosynthetic pathways. (B) A molecular network 
generated from down-regulated genes with the addition of SREBPs to the network. Green lines represent activation and red lines represent 
inhibition. 

doi:1 0.1 371 /journal.pone.01 05674.g006 



GeneGo, the top 10 down-regulated processes belonged primarily 
to cholesterol and lipid biosynthetic pathways (Fig. 6A). As with 
the up-regulated genes, when a molecular network was created 
with the down-regulated genes with the addition of SREBPs to the 
network, a large number of these genes were found to be 



interconnected with SREBP-1 and SREBP-2 (Fig. 6B) indicating 
that this is not an isolated event but are molecular events acting in 
concert to down-regulate cholesterol and lipid biosynthesis. 



PLOS ONE I www.plosone.org 



10 



August 2014 | Volume 9 | Issue 8 | el 05674 



ER Retention of Type MB Procollagen in SIP 




WT (S1P^/') 



S1Pc/fo 





BiP 




Sdf2l1 



G. 



1.E+03 



1.E+02 ■ 



1.E+01 ■ 



O 

■g 1.E+00 ■ 
Q. 

tn 

1.E-01 H 



1.E-02 ■ 



1.E-03 ■ 



1.E-04 '!- 



m 

k 



Ddit3 (5.53 fold) 




Insigl (-2.25 fold) 



m 
I 

o 
00 



m 



o 



m 



m 
+ 
o 



m 



Wild Type 



m 
+ 
o 



m 
+ 
o 
w 



Figure 7. Intact ER stress response in SI P"^*" chondrocytes. (A) RNAfrom the chondroepiphyseal cartilage of El 6.5 WT(lane 1) and S1P'^*° (lane 
2) used in microarray analysis were converted to cDNA and amplified by XBP-1 PCR primers to identify spliced (s) XBP-1 mRNA in a 4.8% 
polyacrylamide gel. (B) XBP-1 was amplified by PCR using XBP-1 PCR primers 3S and 12AS with cDNA derived from El 6.5 WT (lanes 1, 3, 5) and SIP"*" 
(lanes 2, 4, 6) epiphyseal cartilage RNA and the PCR product restriction digested with Psfl which selectively cuts the un-spliced (u) XBP-1 mRNA, and 
the resulting products visualized in a 2% agarose gel. In lanes 1 and 2, the cDNA used are from the same embryos used in (A) and for genome-wide 
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expression profiling. Each lane in lanes 3-6 show analyses fronn RNA pooled from the chondroepiphysis of two different embryos. Thus a total of five 
WTand five SIP"^*" embryos were analyzed. The inverse of the gels are shown in both (A) and (B) to enhance visualization of spliced XBP-1 mRNA. (C- 
F) Expression signaling for two ATF6-driven genes, BiP and Sdf2l1, as seen by in situ hybridization analyses in WT and SI P'^*° cartilage. BiP expression 
is seen in the ulna, carpal, and metacarpal regions in El 6.5 WT (C) and SI P*^*" (D) forelimbs. Sdf2ll expression is seen in the femur of El 5.5 WT (E) and 
j^pcto jpj gg|.. scatter plot generated from quantitative real-time PCR analysis in the murine UPR RT^ Profiler PCR Array system 

comparing the relative expression of 84 genes between WT and SI P"^*" chondrocytes. A log transformation plot is shown in which the relative gene 
expression level of each gene (2"'^'-') in WT is plotted against the corresponding value in SI P'^*° to indicate fold changes (2"'^'^*-'). The black line 
indicates no fold change (fold change of 1 ). The pink lines indicate a fold change of 2 (gene expression threshold). All genes within these two lines are 
considered to be similar in expression to WT. Only Ddit3 or lnsig-1 were significantly differentially expressed among the 84 genes profiled. A total of 
four WT and four SI P*^*" embryos were profiled. 
doi:l 0.1 371/journal.pone.01 05674.g007 



Intact unfolded protein response in siP chondrocytes 

The UPR is modulated by the action of three distinct pathways 
via three distinct ER transmembrane sensors, IREl (inositol- 
requiring transmembrane kinase and endonuclease 1), PERK 
(protein kinase-like ER kinase), and ATF6 [26-28]. Activation of 
IRE 1 results in the endonucleolytic cleavage of XBP- 1 mRNA to 
allow for frame-shift translation and expression of the active XBP- 
1 transcription factor, which increases the expression of UPR- 
responsive genes. Activation of PERK pathway leads to phos- 
phorylation of elongation initiation factor 2 a thereby inactivating 
it, which attenuates translation thus decreasing the protein burden 
of the cell. Furthermore, this diminished translation allows for the 
selective increased translation of ATF4 which alone or in 
combination with XBP-1 allows for expression of UPR-responsive 
genes. Release of ATF6 from the ER by the combined action of 
SIP and site-2 protease (S2P) and its translocation to the nucleus, 
in cooperation with NF-Y or XBP- 1 , up-regulates transcription of 
UPR-responsive genes. 

An indication that the PERK pathway is normal in SlF*" is 
seen by the up-regulation of ATF4 identified in the microarray 
analysis (Fig. 4B, 5B). This is further confirmed by the increased 
expression of ATF4-controlled genes such as DNA damage 
inducible transcript 3 (DditS, also known as C/EBP zeta/ 
CHOP/CHOP 10; Fig. 4B, 5B) [29]. Ddit3 is known to be 
involved in a variety of ER stress response pathways [30,31]. It is 
also known to induce apoptosis in response to ER stress [32]. This 
up-regulation is in consensus with the drastic increase in apoptosis 
seen in SIP^*" cartilage [15]. ATF4 also up-regulates genes 
involved in amino acid metabolism that may protect against 
oxidative stress by promoting glutathione synthesis. The micro- 
array analysis also showed increased expression of two well known 
markers of this response, the amino acid transporter protein solute 
carrier family 7 (Slc7a3) and the amino acid biosynthetic enzyme 
asparagine synthetase (Asns) (Fig. 4B and 5B) [33,34]. Up- 
regulation of se\'eral aminoacyl tRNA synthetases such as seryl 
amino-acyl tRNA synthetase (Sars) and alanyl tRNA synthetase 
(Aars) is also seen (Fig. 4B). Advanced ER stress response is also 
indicated by another target of ATF4-CHOP, namely tribbles 3 
(Trib3/Nipk) (Fig. 4B, 5B), an ER stress inducible gene, which 
also functions by down-regulating its own induction to attenuate 
CHOP and ATF4 mediated transcriptional events [35]. Induction 
of the cell cycle inhibitor p21 (Cdknla) (Fig. 4B, 5B), a negative 
regulator of apoptosis [36] could be the cell's attempt to 
antagonize increased apoptosis. 

To analyze events directed by the IREl pathway we analyzed 
XBP-1 sphcing in S1P'*° chondrocytes (Fig. 7 A, B). Using XBP-1 
primers whose sequences are common to both the spliced (s) and 
un-spliced (u) XBP-1 mRNA [21,25] we checked the ability of the 
primers to generate the smaller, spliced XBP-1 mRNA product in 
PCR reactions using cDNA made from and SIP^*" RNA used 
in the genome-wide expression profiling (Fig. 7A). Both WT and 
SIP"*" cDNA showed the ability to generate the s-XBP-1 mRNA. 



We also confirmed XBP-1 splicing by coupling restriction 
digestion with PCR by using primers 3S and 12AS [22] on cDNA 
derived from the chondroepiphyseal cartilage RNA of more than 
one E16.5 WT and S1F*° embryo (Fig. 7B). Again, both WT and 
S1P^*° showed the ability to splice XBP-1 mRNA suggesting an 
intact IRE-1 directed UPR system in the mutant chondrocytes. 

Based on microarray analysis, we did not observe any down- 
regulation of any known ATF6-responsive genes during ER stress 
response. Genes whose expressions are driven by ATF6 such as 
glucose related protein 78 (Grp78/BiP/Hspa5), heat shock protein 
90bl (Hsp90bl/Grp94), derlin-3, and stromal cell-derived factor 
2-like 1 (SdBU), and are significandy down-regulated during ER 
stress induction in ATF6—/— cell lines [37] were not identified as 
differentially regulated in SIP*" cartilage. To confirm this 
observation we performed in situ hybridization analysis for Sdl211 
and BiP. Significant expression of these genes was detected in 
gjpcto cartilage (Fig. 7D, F). In order to validate the microarray 
results and to get a definitive, global answer about ER stress 
response we performed qPCR analysis by using the murine 
Unfolded Protein Response RT Profiler PCR Array system that 
profiles the expression of 84 genes (Table S2) involved in 
recognizing and responding to misfolded protein accumulation 
in the ER (Fig. 7G). For this assay we used fresh RNA isolated 
from the chondroepiphysis of four WT and four S1P^*° El 6.5 
embr\'os allowing for biological duplicates. Only two significant 
deviations from WT cartilage was observed, namely Ddit3 and 
Insig-1, which were also identified from microarray studies 
(Fig. 4B). The rest of the genes profiled including genes such as 
BiP (HspaS), Herp, and Edeml, whose expressions were down- 
regulated on ER stress in ATF6—/— cell lines, did not show any 
significant change in the mutant. These observations confirmed 
normal UPR responses in S1P^*°. 

Down-regulation of SREBPs-responsive genes in S1P'^'^° 
cartilage 

Figure 6 suggested that genes/pathways regulated by SREBPs 
involved in cholesterol and fatty acid biosynthetic processes are 
primarily down-regulated in S1P^*° cartilage. In order to confirm 
this observation we performed qPCR analyses of a number of 
genes identified as down-regulated (see Table S4). A numbc'r of 
these are known to be direct targets of SIP/SREBP regulatory 
pathway, i.e. genes reported to have a sterol regulatory element 
(SRE) motif in their promoter and requiring SREBPs as 
transcription factors for their synthesis [38] (see Fig. 4B, 6B). As 
with the qPCR-based global UPR analyses above, we used fresh 
RNA pooled from the chondroepiphysis of five El 6.5 WT or five 
SIP'*" embryos. A number of genes identified as down-regulated 
by microarray analysis such as Fgfrlop2, Sparc, and Loricrin did 
not show down-regulation in mutant chondrocytes on qPCR 
analysis (not shown). Along with Prg4 (which is not reported to use 
SREBPs as transcription factors for its expression; see DISCUS- 
SION), qPCR analyses confirmed that genes such as Scdl, 
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Table 2. Quantitative real-time PCR analyses of genes down- 
regulated in S1P'^''° chondrocytes when compared to WT 
littermates. 



Gene 


Wild Type 


SIP'*" 


Scdl 


1.00 


0.27 


Sc4mol 


1.00 


0.41 


Fads2 


1.00 


0.54 


Scd2 


1.00 


0.57 


Stard4 


1.00 


0.58 


Ldir 


1.00 


0.59 


Insig-1 


1.00 


0.63 


Fdps 


1.00 


0.65 


Hmgcs 


1.00 


0.86 


Hmgcr 


1.00 


0.91 



RNA was harvested and pooled from the chondroepiphyseal cartilage of five 
E16.5 SIP'*" or five E16.5 WT embryos. 
doi:l 0.1 371 /journal.pone.Ol 05674.t002 



Sc4mol, Fads2, Scd2, Stard4, Ldlr, Insig-1, and Fdps, which are 
direct targets of SREBPs and involved in the cholesterol/lipid 
biosynthetic pathways, are down-regulated in SlF*" cartUage 
(Table 2) mirroring that seen in the microarray. Especially 
significant among these genes was Scdl which exhibited the 
largest down-regulation. We further validated these observations 
by in situ hybridization analysis demonstrating the down- 
regulation of Scdl, Scd2, Ldlr, Fads2, and Fdps genes at the 
RNA level in the humerus of SlP^*" mice (Fig. 8). Immunohis- 
tochemistry also confirmed the down-regulation of the Scdl 
protein in the distal femoral cartilage in SIP'^'"' mice (Fig. 8). These 
analyses confirmed the significant down-regulation of cholesterol 
and fatty acid biosynthetic pathways in the mutant and suggested 
that this down-regulation may be causal to the abnormal cartilage 
ECM and lack of endochondral bone development. 

To follow this observation further, we studied Col II retention 
and endochondral bone formation in Scdl knockout [Scdl—/—) 
mice. Among genes with a direct molecular link to SREBPs, the 
Scdl gene was maximally down-regulated and thus could be 
considered a good candidate to induce the SIP*" phenotypes. 
Scdl enzyme catalyzes the A9-desaturation of saturated fatty acids 
such as palmitic acid and stearic acid to their corresponding 
monounsaturated fatty acids (MUFA), palmitoleate and oleate, 
respectively. As membranes are composed of lipids, we hypoth- 
esized that retention of matrix proteins in the ER and the resulting 
ER stress could be due to ER membrane lipid changes. Changes in 
MUFA composition of the ER membrane and the resulting 
changes in membrane fluidity could result in Col II retention in 
the ER. To understand whether the phenotype in SIP^*" mice is 
driven by lack of Scdl activity, we analyzed endochondral bone 
formation and Col II deposition in Scdl—/— mice. Scdl — / — 
mice not only suffer from the down-regulation of Scd 1 activity but 
also demonstrated down-regulation in lipid biosynthesis 
[22,39,40], a situation that could be considered analogous to the 
down-regulation of cholesterol and lipid homeostasis seen in 
SIP'*" chondrocytes. Interestingly however, onset of endochondral 
bone formation and Col II deposition in Scdl—/— mice were 
similar to that seen in wild type [Scdl+/+) mice (Fig. 9). 



Discussion 

The disruption of Col II secretion is central to SIP ablations 
[15,41]. In this study we demonstrated that SIP ablation in 
cartilage results in the specific retention of the pro-Col IIB isoform 
in the chondrocyte ER. As the Col2-Cre system initiates deletion 
of SIP very early at E11.5-E12 [42], chondroprogenitor cells 
which synthesize Col IIA also suffer from SIP ablation. However 
this has no effect on Col IIA secretion, a difiFerence probably 
mediated by the presence of exon 2 coded sequences in the N- 
propeptide of Col IIA. Thus, pro-Col IIB is the major cartilage 
protein retained in mutant chondrocytes which prevents the 
formation of a Col II matrix in SIP*" cartilage. This indicates that 
SIP ablation disrupts pro-Col IIB trafficking from the ER and that 
a S IP-derived activity is essential for pro-Col IIB trafficking from 
the ER to the cartilage. The minor retention of COMP is 
presumably due to interactions between small amounts of COMP 
with the large pro-Col IIB aggregates in the ER. This phenom- 
enon is not unique to SlP*° and has been observed in 
pseudoachondroplasia (PSACH) chondrocytes. PSACH chondro- 
cytes express mutant COMP which causes its retention in the ER, 
Interactions between the trapped mutant COMP and Col II 
results in matrix formation inside the ER [43] with Col II 
aggregates surrounded by a COMP matrix. These observations 
indicate that the minor retention of COMP seen in SIP'*" 
cartilage is incidental to pro-Col IIB retention. The dynamic 
association between Col II and COMP in the cartilage is also 
evidenced by the similarity and overlap of their distribution in WT 
cartilage and the absence of a similar COMP organization in 
SlP*" cartilage due to the lack of a pro-Col IIB-derived Col II 
matrix. 

In order to identify the molecular pathways affected by SIP 
ablation that disrupted pro-Col IIB trafficking, we performed 
genome-wide expression analysis in the SIP*" model. Genome- 
wide expression profiling on RNA isolated from the chondroepi- 
physeal cartilage of El 6.5 SlP*° and WT littermates followed by 
2-way ANOVA and gene ontology enrichment analysis demon- 
strated that only the SREBPs-regulated pathways are affected in 
SlP*" chondrocytes. The Prg4 (lubricin) gene which has no 
reported association with SIP or SREBPs or any SRE motif in its 
promoter was also found to be down-regulated. However, 
knockout mouse models for Prg4 [44] do not show any defects 
in endochondral bone formation or Col II trafficking to the ECM. 
Therefore, it seems possible that Prg4 down-regulation in SIP'*" 
cartilage is a secondary effect due to poor cartilage matrix 
development in these mice or due to abnormal chondrocyte 
differentiation. 

As chondrocytes secrete large amounts of various cartilage 
matrix proteins to assemble the specialized cartilage ECM, a 
functional UPR would be indispensable to this function. Thus we 
had hypothesized that SIP ablation would prevent an effective 
UPR response to ER stress which would cause intracellular Col II 
accumulation [15]. However, UPR response pathways are 
unaffected in SlP*° chondrocytes. The activity of ATF6 is 
unchanged on SIP ablation and there is a lack of differential 
expression of ATF6-responsive genes such as BiP and Sdf211. This 
can be explained by reasoning that the requirement of SIP for 
ATF6 processing is not absolute, which is also noted in other 
systems. ER stress response was reported to be normal in zebrafish 
gonzo phenotype where expression of BiP was shown as unaltered 
in absence of SIP [6]. The S IP-lacking SRD-12B cells demon- 
strated only a partial ability to prevent ATF6 processing to its 
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Figure 8. Expression analyses of some key genes involved in cholesterol and fatty acid biosynthesis in SIP*^*" cartilage by in situ 
hybridization. Notice the drastic reduction in the expression of Scdl, Scd2, Ldir, Fads2, and Fdps genes at the RNA level, shown by in situ 
hybridization in SI P^*° humerus and by immunohistochemistry (IHC) for Scdl protein in SI 9^''° femur. All analyses are shown at El 5.5 except for Fdps 
which is shown at El 4.5. Bars: (in situ hybridization analyses): 500 |a.m; (IHC): 200 |im. 
doi:1 0.1 371 /journal.pone.01 05674.g008 



active form on induction of ER stress; only S2P ablation allowed 
for a complete lack of ATF6 processing [45]. Thus, even moderate 
amounts of active ATF6 in concert with other nuclear factors such 
as NF-Y [30] and/or XBP-1 [46] would have sufficient 
transcriptional induction properties to induce UPR. Thus, in this 
study it was not possible to address the importance of UPR to 
cartilage ECM development. However, viable and fertile ATF6a 



knockout mice were reported [46] suggesting that ATF6-directed 
ER stress responses are not necessary for normal Col II deposition 
or endochondral bone development in mice. Notably, these studies 
indicate that the apoptosis seen in SIP^*" cartilage is due to the 
trapped Col II in the ER and induction of UPR, rather than an 
inabihty to respond to ER stress. 
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Figure 9. Normal endochondral bone development and Col II 
deposition in Scdl—I— mice. (A, B) Sections of El 5.5 humerus were 
stained with Safranin O, Fast green, and hematoxylin showing normal 
onset of endochondral bone development in Scdl—/— mice. (C-F) 
Double-labeled immunofluorescence analyzing Col II deposition in 
Scdl—/— mice using IIP and IIA antibodies. Colors represent antibody 
localizations as follows: green, Col IIA (IIA antibody); red, Col II THD (HP 
antibody); yellow, colocalization of both antibodies; blue, DAPI-stained 
nuclei. Panels C and D show the matrix around early immature 
chondrocytes in the resting zone; panels E and F show the matrix 
around mature columnar chondrocytes in the proliferative zone. Bar: (A, 
B): 500 iim; (C-P): 50 ^m. 
doi:1 0.1 371/journal.pone.01 05674.g009 

Other ER stress transducers such as OASIS and its related 
family members BBF2H7 and CREBH also appear to be 
unafiected. OASIS has a specific expression pattern and is 
primarily induced in bone tissues and the central nervous system, 
but no expression in chondrocytes has been reported [47]; likewise 
CREBH is liver-specific [5]. Thus OASIS and CREBH and their 
dependent pathways are not relevant in S1F^*°. But BBF2H7 
differs from OASIS in its strong expression in the proliferating 
zone of the cartilage [48]; it has also been reported to need SIP for 
its processing [49] and would be expected to be inactive in SIP'*" 
chondrocytes. Bbf2h7-/- mice show severe chondrodysplasia and 
the retention of Col II in the ER similar to that seen in SlK*° mice 
[48]. BBF2H7 is required for the expression of Sec23a protein, a 
component of COPII vesicles which transports secretory proteins 
from ER to the matrix. However, Sec23a is unaffected in SIP'*" 
chondrocytes as it was not identified as diflFerentiaUy down- 
regulated in the microarray analysis and was further validated by 
qPCR analysis (Figure S2); nor were its protein levels affected (not 
shown). Therefore, in S1P^*° chondrocytes, ER stress transducers 
and their activities have remained largely unaffected. 

However the requirement for SIP is absolute for the SREBPs. 
SRD-12B (SIP-/-) ceUs, which only partially prevent ATF6 
processing, are able to impose a complete block of SREBPs 



processing [4,5, .50]. This would explain the specific down- 
regulation of SREBP-driven pathways in SIP'*" chondrocytes. 
However, this down-regulation does not appear to be responsible 
for the mutant phenotypes seen on SIP ablation in chondrocytes. 
We had surmised that potential changes in ER membrane 
composition and fluidity accompanying the down-regulation of 
sterol and lipid pathways may be responsible for Col II retention 
seen in the ER. If this was true then these ER membrane changes 
would also be present in the chondrocyte ER of Scdl —/— mice 
that would lead to abnormal cartilage and endochondral bone 
development. However, we observed no intracellular Col II 
retention or lack of endochondral bone development in Scdl —/— 
mice. It is also striking that mice lacking Ldlr [51,52] or Fads2 
(Delta-6 desaturase) [53], genes that are also down-regulated in 
SIP'*" chondrocytes, show no disruption in endochondral bone 
development as that seen in SIP'*". Furthermore, neither are these 
mutant phenotypes seen in mouse models with compound or 
multiple disruptions in lipid homeostasis such as Ldlr—/—;Lcat—/ 
— mice [51]. 

These observations may be explained by the possibility to 
separate cartilage developmental pathways from lipid pathways. In 
the zebrafish gonzo phenotype, SIP ablation results in both 
cartilage and lipid phenotypes. But ablation of SREBP cleavage- 
activating protein (SCAP) results only in lipid phenotypes but no 
cartilage phenotype [6]. This indicates that lipid and cartilage 
defects are caused by different mechanisms. However, it has been 
reported that the inhibition of cholesterol biosynthesis through the 
administration of a chemical inhibitor AY 9944 in mice suppressed 
longitudinal bone growth via suppression of chondrocyte prolif- 
eration and hypertrophy [54]. These contrasting observations 
make it difficult to discern the exact mechanistic input of sterol and 
lipid homeostasis to cartilage development. It is possible that 
cholesterol/lipid down-regulation and ER entrapment of Col II 
are independent phenotypes from two separate functions of SIP, 
but those functions relating to Col II trafficking remain to be 
identified. 

Figure 10 summarizes our current findings. In this study we 
investigated all current known functions of S 1 P: induction of ER 
stress response, induction of cholesterol and fatty acid biosynthesis, 
processing of BBF2H7 to effect Sec23a expression, and the 
processing of the ot/ P subunit precursor of N-acetylglucosamine- 1 - 
phosphotransferase complex (GNPTAB). Recently it was demon- 
strated that GNPTAB, involved in the addition of mannose 6- 
phosphate residues on lysosomal enzymes, is a direct SIP substrate 
[55]. In our preliminary experiments, there is at least equivalent 
GNPTAB ot/ fi precursor cleavage activity in S 1 P'^*" chondrocytes 
as in chondrocytes from its wild type littermates (not shown). This 
indicates that this activity is not lost in SIP^*" chondrocytes and 
therefore not causal to the mutant phenotypes. Besides, mucolip- 
idosis II {Gnpta¥-^"^^""^-) mice (which lack GNPTAB) do not 
exhibit lack of endochondral bone formation or Col II entrapment 
in the ER [55]. Our studies thus indicate that activities associated 
with current known functions of S 1 P cannot provide for a direct 
molecular explanation of SIP'^*" phenotypes, especially the 
entrapment of pro-Col IIB in the ER. We therefore propose that 
S 1 P has additional, novel substrates in chondrocytes (schematized 
in Fig. 10), likely non-transcription factor-related, that is vital to 
normal cartilage development; at least one of these substrates is 
specific and indispensable for pro-Col IIB processing and/or 
trafficking from the ER. The identification of this novel SIP 
substrate will not only allow for a molecular understanding of 
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Figure 1 0. A schematic summarizing the results from this study and possible additional roles for SI P in cartilage development. S1 P 

is historically known for processing the transcription factors ATF6 and SREBPs. Two relatively recently discovered activities include the processing of 
BBF2H7 and the ot/p subunit precursor of N-acetylglucosamine-1 -phosphotransferase complex (GNPTAB). Of these four known activities, only SREBPs 
activities relating to cholesterol and fatty acid homeostasis are down-regulated in SIP"^*" chondrocytes, with maximal down-regulation of Scdi. 
However this does not appear to be responsible for the SI 9'^''° mutant phenotypes as pro-Col MB retention in chondrocytes or lack of endochondral 
bone development is not seen in Scdl ~/~- mice. The retention of pro-Col MB in the ER could be due to changes in ER membrane lipid composition, 
though unlikely as pro-Col MB retention is not seen in Scdi —/— mice. Processing of additional, yet unidentified, SI P substrates in chondrocytes likely 
modulate pro-Col MB trafficking from the ER. Lack of endochondral bone formation is probably due to the abnormal cartilage matrix devoid of Col II, 
or due to unidentified SI P-regulated processes. In the figure, blue arrows indicate activities that are normal in S1P'^''° chondrocytes; red arrows 
indicate activities that are down-regulated or abnormal in SIP^*". Dashed arrows with a question mark indicate the possibility of these mechanisms 
directed by novel SIP substrates, but nothing more is known as yet. 
doi:10.1371/journal.pone.0105674.g010 



SIP*" phenotypes, but will also allow for an unraveling of the 
molecular mechanism involved in pro-Col IIB trafficking, a 
phenomenon that still remains poorly understood. Our studies 
thus indicate an increased breadth of S 1 P functions in chondro- 
cytes mandatory for normal cartilage development. 

Supporting Information 

Figure SI (A) A diagrammatic representation of the W T 
embryonic growth plate. All immunohistochemical analyses 
are shown for the columnar proliferative regions above the 
hypertrophic cells and for a corresponding region in S1P^*°. The 
chondroosseous junction is the approximate region where the tips 
of long bones were broken off and harvested for cartilage analysis 
by microarray . S 1 P^*° lacks endochondral bone and the diaphysis 
of the bone is cartilaginous. (B) A diagrammatic representation of 
the cartilage tissue harvested for microarray analysis. For analysis 
of WT and S 1 P^*° cartilage by microarray, the tips of the long 
bones (shown demarcated by the dashed line and denoted as 
epiphyseal cartilage) were broken off and stored in Trizol. When 
entire long bones were harvested for analysis, the bones were 
harvested intact and the epiphyseal cartilage portions were not 



separated from the bone. (C) A schematic representation of the 
differences between WT (rectangles) or S1P^*° (triangles) RNA 
using a principal component analysis scatter plot based on 
microarray hybridization signals. RNA was analyzed from 
cartUage (green), entire long bones (red), and calvaria (blue). 
Differential gene expression between WT and SIF*" calvaria 
(blue) is minimal (located very close to each other). Differential 
gene expression between WT and SIP^*" cartilage (green) is 
maximum (located farthest from each other), while those between 
WT and S1P^*° long bones (red) is intermediate to cartUage and 
calvaria. 
(EPS) 

Figure S2 Expression levels for Sec23a are unaffected in 
SIP"*" mice. Each WT or SIP'*" (CKO) analysis used RNA 
pooled from the chondroepiphyseal cartilage of two embryos. 
Thus a total of four WT and four S 1 P'^*" embryos were analyzed. 



Q;PCR analysis was done using the 2 ' 

values were set to one. 

(EPS) 



method where the WT 



Table SI SYBR Green qPCR primer sets for the murine 
genes listed were taken from the MGH/Harvard 
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Medical School primer bank (http://pga.mgli.liarvard. 
edu/primerbank/) . 

(DOCX) 

Table S2 A complete list of 84 genes that were profiled 
by qPCR using the murine Unfolded Protein Response 
RT^ Profiler PGR Array system, shown grouped accord- 
ing to their functions. 

(DOCX) 

Table S3 A partial list of genes significantly up- 
regulated in SIP"*" chondrocytes identified from micro- 
array analysis. Shown are top 20 genes that are differentially 
up-regulated when compared to WT, and some selected genes. 
RNAs with no corresponding gene names are not included in this 
list. 

(DOCX) 

Table S4 A partial list of genes significantly down- 
regulated in S1P^*° chondrocytes identified fi-om micro- 
array analysis. Shown are top 20 genes that are differentially 
down-regulated when compared to WT, and some selected genes 
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